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Abstract 

The  viscous,  three-dimensional,  flowfields  of  UH60 
and  BERP  rotors  are  calculated  for  lifting  hover  con¬ 
figurations  using  a  Navier-Stokes  computational  fluid 
dynamics  method  with  a  view  to  understand  the  im¬ 
portance  of  planform  effects  on  the  airloads.  In  this 
method,  the  induced  effects  of  the  wake,  including  the 
interaction  of  tip  vortices  with  successive  blades,  are 
captured  as  a  part  of  the  overall  flowfield  solution  with¬ 
out  prescribing  any  wake  models.  Numerical  results  in 
the  form  of  surface  pressures,  hover  performance  pa¬ 
rameters,  surface  skin  friction  and  tip  vortex  patterns, 
and  vortex  wake  trajectory  are  presented  at  two  thrust 
conditions  for  UH60  and  BERP  rotors.  Comparison  of 
results  for  the  UH60  model  rotor  show  good  agreement 
with  experiments  at  moderate  thrust  conditions.  Com¬ 
parison  of  results  with  equivalent  rectangular  UH60 
blade  and  BERP  blade  indicates  that  the  BERP  blade, 
with  an  unconventional  planform,  gives  more  thrust  at 
the  cost  of  more  power  and  a  reduced  Figure  of  Merit. 
The  high  thrust  conditions  considered  produce  severe 
shock-induced  flow  separation  for  UH60  blade,  while 
the  BERP  blade  develops  more  thrust  and  minimal  sep¬ 
aration.  The  BERP  blade  produces  a  tighter  tip  vortex 
structure  compared  with  the  UH60  blade.  These  results 
and  the  discussion  presented  bring  out  the  similarities 
and  differences  between  the  two  rotors. 


Nomenclature 


doo  — 

speed  of  sound,  characteristic  velocity 
scale 

6  = 

number  of  blades 

c(r)  = 

blade  chord  at  radial  location  r 
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<^eq 

= 

equivalent  blade  chord,  1/R  f  c{r)dr 

Cttu 

=r 

thrust- weighted  equivalent  chord, 
f  c(r)r^dr/ f  r~dr 

Cp 

pressure  coefficient  bsised  on  local 
dynamic  pressure 

Cq 

= 

rotor  torque  coefficient,  Q/(2q(R)5rR^) 

Ct 

sectional  thrust  coefficient  based  on 
local  chord,  dt/dr/(c(r)q(R)) 

sectional  thrust  coefficient  based  on 

Ce„  cff/dr/(ce,q(R)) 

Ctavt 

= 

area  under  Ct  vs  r  curve 

Ctavt 

= 

area  under  vs  r  curve 

Ct 

= 

rotor  thrust  coeffient,  T/(2q(R)n-R^) 

e 

= 

energy  per  unit  volume 

E,F,G 

= 

flux  vectors 

FM 

Figure  of  Merit,  /Cq 

Mtip 

= 

blade  tip  Mach  number 

q{R) 

= 

rotational  dyneunic  pressure. 

1/2  p{n  R)2 

Q 

= 

flow-field  vector 

Q 

= 

rotor  torque 

r 

rotor  radial  coordinate 

R 

= 

blade  radius 

Re 

= 

Reynolds  number  based  on  Mup  and 

/» 

S 

= 

Ceq 

viscous  flux  vector 

T 

rotor  thrust 

U,  V,  w 

= 

velocity  components  in  physical  space 

^loc 

local  blade  speed 

x,y,z,t 

■  = 

physical  space  coordinates 

a 

z= 

angle  of  attack,  degrees 

r 

= 

sectional  circulation 

f 

= 

normalised  blade  sectional  circulation, 
T/QR? 

0. 

= 

collective  pitch 
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p 

<T 


(Tx 


n  = 


density 

rotor  geometric  solidity, 
bf  c(r)rdr/(7rR^) 
rotor  thrust-weighted  solidity, 
b  /  ctwrdr/{irR?) 

generalized  curvilinear  coordinates 
rotor  angular  velocity 

Introduction 


The  accurate  prediction  of  the  flowfield  of  a  he¬ 
licopter  rotor  continues  to  be  one  of  the  most  com¬ 
plex  and  challenging  problems  in  applied  aerodynam¬ 
ics.  This  is  true  despite  the  availability  of  the  present 
day  supercomputers  of  Cray-2  class  and  improved  nu¬ 
merical  algorithms.  The  complexity  of  the  flow  stems 
from  several  peculiar  features  that  are  unique  to  a  heli¬ 
copter  rotor;  the  dominant  of  these  is  the  vortical  wake, 
which  strongly  influences  the  operating  characteristics. 
Accurate  prediction  of  this  vortical  wake  is  probably 
the  most  important,  most  studied  and  the  most  diffi¬ 
cult  aspect  of  helicopter  flowfield. 


Current  methods  of  analysis  of  wakes  range  in  com¬ 
plexity  from  relatively  simple  momentum-theory  ap¬ 
plications  to  sophisticated  free-wake  methods.  In  be¬ 
tween  these  extremes,  there  exists  a  variety  of  so-called 
prescribed- wake  models^  Although  such  methods  are 
often  used  in  industry  for  predicting  airloads,  they  do 
suffer  from  the  limitation  that  the  empirical  determi¬ 
nation  of  wake  shape  ignores  some  of  the  important 
flowfield  details  such  as  the  mutual  interaction  of  in¬ 
dividual  vortex  elements.  The  fact  that  most  of  these 
methods  have  to  be  coupled  with  some  wake  model,  for 
realistic  estimates  of  airloads,  restrict  their  application 
to  blade  shapes  and  planforms  that  are  geometrically 
simple.  Besides,  further  modeling  is  required  to  simu¬ 
late  the  formation  of  concentrated  tip  vortex  in  the  tip 
region. 


The  recent  application  of  Computational  Fluid  Dy¬ 
namics  (CFD)  methods  to  predict  these  complex  flows 
has  opened  a  new  avenue  to  address  these  problems 
more  rigorously.  Methods  based  on  the  solutions  of 
full-potential,  the  Euler  and  the  Navier-Stokes  equa¬ 
tions  abound  in  the  literature.  Although  most  of  these 
methods  still  use  some  sort  of  wake  model  to  bring-in 
the  influence  of  vortical  wake,  there  are  methods  that 
capture  the  effects  of  wake  fairly  accurately  as  a  part  of 
the  overall  flowfield.  In  this  regard,  the  recent  methods 
based  on  the  solutions  of  full-potential^,  the  Euler^’^ 
and  the  Navier-Stokes®  equations  are  noteworthy. 


In  the  present  study,  the  Navier-Stokes  methodol¬ 
ogy  of  Srinivasan  et  al.®  is  used  to  calculate  the  flow- 
fields  of  realistic  helicopter  rotor  configurations  of  the 
UH60®  and  BERP^  rotors  for  lifting  hover  configura¬ 
tions,  with  a  view  toward  assessing  the  importance 
of  planform  effects  on  the  airloads.  Calculations  are 
also  performed  for  equivalent  rectangular  blade  of  the 
UH60  rotor  to  further  understand  the  role  of  planform 
in  determining  the  airloads  and  tip  vortex  wake  geom¬ 
etry.  The  numerical  code  of  Ref.  5,  called  TURNS 
(Transonic  Unsteady  Rotor  Navier-Stokes),  will  be  uti¬ 
lized  for  this  purpose.  This  methodology  has  been  val¬ 
idated  with  experimental  data  for  aerodynamics  and 
aeroacoustics  of  high-speed  impulsive  noise  for  rotors 
in  hover  and  forward  flight*.  Currently  the  vortical 
structure  calculated  with  this  method  is  smeared  due 
to  grid  coarseness.  However,  the  nearfield  wake  tra¬ 
jectory  is  captured  fairly  accurately.  The  use  of  the 
Navier-Stokes  method  enhances  (a)  the  tip  flow  sim¬ 
ulation  which  involves  resolving  the  three-dimensional 
separated  flow  and  concentrated  tip  vortex,  and  (b)  the 
accurate  simulation  of  strong  viscous-inviscid  interac¬ 
tion  involving  shock-induced  separation  at  high  blade 
tip  speeds. 

Covering  Equations  and  Solution  Methodology 

The  governing  partial  differential  equations  are 
the  unsteady,  thin-layer  Navier-Stokes  equations.  For 
generality,  these  equations  are  written  in  strong  conser¬ 
vation-law  form  in  gereralized  body-conforming  curvi¬ 
linear  space  (^,  77,  (^,  r)  and  are  given  by® 


drQ  +  d^E  +  d^F^dQG  =  Re-^d(,S  (1) 

where 

Q  =  \p,pu,pv,pw,ef,Q  =  Q/J  (2) 


is  the  flowfield  vector  being  solved  for  and  E,  F,  G  are 
convective  flux  vectors  and  S  is  the  viscous  flux  vector 
being  considered  in  the  thin  layer  limit;  J  is  the  Jaco¬ 
bian  of  transformation.  Re  is  the  Reynolds  number  and 
the  sign  ^indicates  that  the  quantity  is  normalized  by 
the  Jacobian.  The  primitive  variables  of  Eq.  (1)  are 
the  density  p,  the  three  mass  fluxes  pu,  pv  and  pw  in 
the  three  coordinate  directions  i,  y  and  z,  respectively 
and  the  energy  per  unit  volume  e.  All  these  quantities 
are  nondimensionlized  by  the  appropriate  free  stream 
reference  quantities.  The  chord  of  the  rotor  blade  and 
the  free  stream  sound  speed  are  assumed  to  be  the  char¬ 
acteristic  length  and  velocity  scales.  The  equation  set 
given  by  Eq.  (1)  together  with  the  equation  of  state  for 
a  perfect  gas  describe  the  entire  flowfield. 
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The  finite-difference  implicit  numerical  scheme  for 
the  solution  of  Eq.  (1)  is  described  in  Ref.  5.  Briefly, 
the  numerical  algorithm  uses  a  Roe  upwind-biased 
scheme  for  all  three  coordinate  directions  with  higher- 
order  MUSCL-type  limiting  for  the  right-hand  side  to 
model  shocks  and  for  propagating  acoustics  accurately. 
The  use  of  Roe’s  upwinding  eliminates  the  addition  of 
explicit  dissipation  and  thus  produces  a  less  dissipa¬ 
tive  numerical  scheme.  The  LU-SGS  (lower-upper  sym¬ 
metric  Gauss-Seidel)  implicit  operator  used  for  the  left- 
hand  side  makes  the  numerical  scheme  computationally 
efficient  eis  well  as  robust.  The  option  of  Newton  sub¬ 
iterations  at  each  time-step  reduces  the  linearization 
and  factorization  errors  as  well  as  makes  the  numerical 
scheme  second-order  accurate  in  time.  This  is  partic- 
ulary  useful  for  unsteady  calculations  (forward  flight). 
Currently  the  scheme  is  third-order  accurate  in  space 
and  first-  or  second-order  accurate  in  time. 

To  take  advantage  of  the  quasi-steady  nature  of 
the  hovering  rotor  flowfield,  the  governing  equations 
are  solved  in  the  blade-fixed  coordinate  system.  This 
necessitates  addition  of  source  terms  to  account  for  the 
centrifugal  force  of  rotation  of  the  blades.  The  coordi¬ 
nate  transformation  space  metrics  are  evaluated  using 
a  finite  volume  formulation.  Although  such  a  method 
captures  free-stream  accurately,  formation  of  time  met¬ 
rics  rigorously  with  this  scheme  is  computationally  ex¬ 
pensive.  So,  a  finite-difference  formulation  is  adopted 
for  unsteady  metrics.  Further,  a  locally-varying  time- 
step  is  used  in  the  integration  procedure  to  accelerate 
convergence^®  of  numerical  scheme.  With  these  fea¬ 
tures,  the  code  is  robust,  fairly  efficient,  and  runs  at 
140  MFLOPS  on  a  Cray  Y-MP  supercomputer. 

Body  conforming  finite-difference  grids  are  used  for 
the  long  slender  blades  of  the  UH60  and  Lynx  BERP 
rotors.  These  blades,  shown  in  planform  in  Fig.  1, 
have  aspect  ratios  of  15.51  and  12.46  respectively.  The 
b  are  highly  twisted  and  have  a  twist  distribution 
alou^  their  lengths  as  shown  in  Fig.  2.  The  charac¬ 
teristics  of  UH60  blade  is  that  it  has  a  hook-like  twist 
distribution  in  the  tip  region  with  a  swept-back  plan- 
form  in  this  area.  The  chord  of  the  blade  is  constant 
along  its  entire  length  with  a  rectangular  planform  ex¬ 
cept  in  the  swept-back  part  of  the  blade.  The  BERP 
blade  on  the  other  hand,  has  a  large  paddle-type  plan- 
form  in  the  outer  part  of  the  blade  near  the  tip.  The 
twist  change  is  linear  like  the  UH60  blade  up  until  the 
paddle  part  where  it  decreases  to  a  negative  constant 
value. 

Because  of  the  axial  symmetry  of  the  hovering  rotor 
flowfield,  a  C-H  grid  topology  is  chosen  with  C-type  in 


the  chordwise  direction  and  H-type  in  the  spanwise  di¬ 
rection.  The  H-type  spanwise  topology  creates  a  bevel 
tip  for  the  square-tipped  UH60  blades.  To  adapt  to 
the  cylindrical  flowfield  of  the  rotor,  three-dimensional 
grids  are  constructed  for  these  blades  by  stacking  and 
bending  two-dimensional  C-H  grids  so  that  they  lie  on 
the  blade  along  the  chord.  Typical  fine  grids  have 
nearly  940,000  points  with  217  points  in  the  chordwise 
direction  and  71  and  61  points  in  the  spanwise  and 
normal  directions,  respectively.  These  grids  have  clus¬ 
tering  in  the  leading  and  trailing  edge  regions  as  well  as 
in  the  tip  region.  The  grid  lines  are  nearly  orthogonal 
at  the  surface  with  a  spacing  of  0.00004  chord  in  the 
normal  direction.  The  grid  boundary  is  located  at  two 
rotor  diameters  at  the  top  and  bottom  boundaries  and 
one  radius  beyond  the  tip  in  the  spanwise  direction. 
Coarse  grids  have  been  constructed  by  eliminating  ev¬ 
ery  other  point  in  all  three-directions  and  typically  have 
120,000  points.  Figure  3  shows  the  views  these  grids 
(for  single  blades)  in  the  plane  of  the  blades.  Figure  4 
shows  a  view  of  a  grid  for  single  blade  constructed  for 
BERP  blade.  This  typical  grid  shows  the  grid  bound¬ 
aries  and  the  details  of  grid  clustering  near  the  grid 
surface  and  in  the  tip  region. 

The  boundary  conditions  are  applied  explicitly.  At 
the  surface  a  noslip  condition  is  used,  but  the  surface 
will  have  finite  velocity  due  to  the  blade  rotation  intro¬ 
duced  through  the  time  metrics.  Characteristic-type 
boundary  conditions  are  used  at  the  farfield;  this  point 
is  discussed  in  more  detail  later  in  the  paper.  The  sur¬ 
face  pressure  is  determined  by  solving  the  normal  mo¬ 
mentum  equation  and  the  density  is  evaluated  by  the 
adiabatic  wail  condition.  At  the  inboard-plane  bound¬ 
ary,  near  the  eixis  of  rotation  of  the  blade,  the  flow 
quantities  are  extrapolated  from  the  grid  interior.  For 
a  hovering  rotor  with  identical  blades,  the  flow  associ¬ 
ated  with  each  blade  is  assumed  to  be  the  same.  To 
take  advantage  of  this  situation,  only  the  flowfield  of 
a  single  blade  is  calculated  by  bringing  information  of 
neighboring  blades  through  periodic  boundaries'*  shown 
in  Fig.  4.  The  flow  of  the  entire  rotor  is  then  assem¬ 
bled  through  a  post-processing  procedure  by  combining 
the  flowfields  of  individual  blades.  Finally,  the  calcula¬ 
tions  reported  here  assume  fully  turbulent  flow  for  the 
blades  and  laminar  flow  in  the  wake.  The  simple  al¬ 
gebraic  turbulence  model  of  Baldwin  and  Lomax**  is 
used  to  estimate  the  eddy  viscosity. 

Resu*»s 

The  present  numerical  scheme  with  its  wake- 
capture  feature  has  performed  better  in  resolving  the 
tip  flow,  shocks  and  overall  airloads  than  many  Navier- 
Stokes  methods  that  have  used  wake  models  or  some 
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variance  of  such  models^^’^^.  In  contrast  to  the  the 
model  two-bladed  hovering  rotor^^  considered  for  calcu¬ 
lations  with  this  solution  method  in  a  previous  study® , 
the  present  study  considers  four-bladed  model  UH60  ro¬ 
tor  and  actual  BERP  rotor.  The  blades  of  these  rotors, 
as  previously  noted,  have  large  aspect  ratio  and  are 
highly  twisted  along  their  lengths  (see  Fig.  2).  In  addi¬ 
tion  they  have  non-rectangular  planforms.  As  before®, 
the  calculations  are  performed  on  isolated  rotor  blades 
and  composite  solutions  are  assembled  by  combining 
the  flowfields  of  all  four  blades.  In  the  calculation  pro¬ 
cedure,  the  evolution  of  the  flowfield  to  a  quasi-steady 
state  is  monitored  as  the  blade  is  set  in  motion  from 
a  quiescent  flow  condition.  Typical  numerical  results 
presented  consists  of  surface  pressures,  sectional  thrust 
distributions,  surface  oil  flow  patterns,  tip  vortex  de¬ 
tails,  and  wake  trajectory. 

UH60  Blackhawk  Rotor 

Figure  5  show  typical  plots  of  surface  pressure  dis¬ 
tributions  for  the  model  UH60  blade  at  several  radial 
stations.  Although  model  hover  tests  are  performed  for 
this  configuration  and  some  plots  of  the  experimental 
data  are  published®  for  this  rotor,  much  of  the  experi¬ 
mental  information  remains  restricted.  In  view  of  this, 
a  representative  tip  speed  {Mup  =  0.628)  and  thrust 
condition  {Ct/<^  =  0.085)  were  selected  and  several 
collective  pitch  settings  (0c)  were  tried  out  to  match 
with  experimental  data.  Results  for  6c  =  9®  showed 
the  closest  match  with  experimental  surface  pressures. 
Figure  5  shows  the  comparison  of  these  results  with 
experiments.  The  numerical  results  appear  to  be  in 
general  agreement  with  the  data  at  all  span  stations. 
The  calculated  Ct/o"  =  0.084  for  this  case  is  also  in 
good  agreement  with  the  experimental  value  of  0.085. 

In  order  to  assess  the  importance  of  planform  ef¬ 
fects,  flowfield  calculations  were  made  on  an  equiva¬ 
lent  rectangular  blade  of  identical  twist  distribution  as 
the  UH60  blade  (shown  as  Twist  1  in  Fig.  2)  for  the 
above  flow  conditions  of  Mnp  =  0.628,  and  dc  =  9® 
and  Re  =  2.75  x  10®.  The  results  are  compared  with 
the  UH60  blade  in  Fig.  6  in  a  plot  of  sectional  thrust 
coefficient  (Ct)  variation  along  the  blade  radius  (r/R). 
In  order  to  compare  these  results  with  experimental 
data®,  the  sectional  thrust  coefficient  has  been  normal¬ 
ized  by  Ctave  which  is  the  area  under  the  curve  of  Ct  vs 
r/R.  Since  the  equivalent  rectangular  blade  has  identi¬ 
cal  twist  ditribution  as  the  UH60  blade  along  its  radius 
and  has  an  unswept  planform  in  the  tip  region,  the 
differences  seen  with  UH60  blade  in  Fig.  6  are  only 
in  this  region.  Although  both  these  blades  have  bev¬ 
elled  tips  (because  of  H-grid  topology),  the  UH60  blade 
produces  a  relatively  smaller  separation  in  this  region 


than  the  equivalent  rectangular  blade  (with  unswept 
tip).  The  tip  vortex  location  for  the  UH60  biade  is 
slightly  inboard  of  that  seen  for  the  equivalent  rectan¬ 
gular  blade.  This  observation  is  in  general  agreement 
with  some  previous  studies  of  swept  and  unswept  plan- 
forms  with  bevel-tips^®.  In  terms  of  performance,  the 
equivalent  rectangular  blade  produces  about  8%  less 
thrust  than  the  UH60  blade,  and  the  Figure  of  Merit  is 
lower.  Table  1  details  the  performance  parameters  for 
these  rotors. 

Figure  7  shows  the  sectional  thrust  distribution 
along  the  radius  of  UH60  blade  for  two  thrust  con¬ 
ditions  corresponding  to  6c  =  9®  and  13°  at  Mup  = 
0.628  and  Re  =  2.75  x  10®.  The  high  thrust  condition 
of  6c  =  13®  produces  a  bound  circulation  distribution 
that  gives  67%  more  thrust  than  the  6c  =  9®  case.  The 
power  requirement  increases  by  a  factor  of  115%.  It  also 
produces  a  much  stronger  tip  vortex  which  is  located 
more  inboard  than  the  6e  =  9®  case.  The  flow  topology 
is  very  complicated  on  the  biade  for  this  high  thrust 
condition  with  shock-induced  boundary  layer  separa¬ 
tion  occuring  over  a  large  part  of  the  blade,  as  shown 
by  the  surface  particle  tracers  (skin  friction  lines)  in 
Fig.  8.  In  contrast,  as  seen  in  the  close-up  view  of  the 
tip  in  Fig.  9,  the  low  thrust  case  {6c  —  9°)  does  not  pro¬ 
duce  shock-induced  flow  separation.  Both  thrust  con¬ 
ditions  produce  flow  separation  on  the  inboard  parts  of 
the  blade  as  seen  in  Fig.8,  but  the  spanwise  extent  is,  of 
course,  greater  at  0c  =  13®.  Figure  10  shows  the  extent 
of  supercritical  flow  on  the  upper  surface  of  the  blade 
in  the  Mach  number  contours  picture  for  the  6c  =  13® 
case.  Also  shown  in  this  figure  are  the  Mach  contours 
in  a  plane  at  r/R  =  0.94.  This  shows  the  location  of 
shock  wave  in  the  leading  edge  region,  and  the  extent 
of  supersonic  pocket. 

Two  representative  performance  curves  of  the 
model  UH60  rotor  as  calculated  here  for  Mup  =  0.628 
are  shown  in  Fig.  11  and  Fig.  12  along  with  UTRC 
data®.  The  details  of  calculated  results  are  shown  in 
Table  1.  In  presenting  these  results,  two  observations 
have  been  made.  First,  the  moderate  thrust  condition 
results  are  in  good  agreement  with  experimental  data. 
Second,  the  high  thrust  condition  considered  in  this 
study  is  apparently  beyond  the  model  test  range.  The 
calculated  Figure  of  Merit  is  approximately  the  same  at 
0c  =  13®  as  at  0c  =  9®  indicating  that  the  rotor  hover 
performance  seems  to  stay  relatively  uniform  for  this 
rotor,  atleast  over  this  range  of  0c. 

In  addition  to  the  planform  change  study  for  this 
blade,  a  fictitious  UH60  blade  having  the  same  plan- 
form,  but  with  linear  twist  in  the  tip  region  (Twist2  in 
Fig.  2),  instead  of  a  hook-type  of  twist  (Twist  1),  was 
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considered  to  examine  the  importance  of  twist  in  this 
region  on  airloads.  The  sectional  thrust  distribution 
for  this  blade  is  shown  in  Fig.  6  along  with  results  for 
the  standard  UH60  blade  at  the  same  test  conditions 
as  of  Fig.  5.  As  seen  in  this  figure,  the  blade  gener¬ 
ates  lower  thrust  inboard  (perhaps  due  to  the  slightly 
lower  values  of  twist  in  this  region),  higher  thrust  at 
r/R  ~  0.8  and  lower  thrust  at  the  tip,  compared  to 
the  UH60  blade.  Overall,  this  blade  produces  15%  less 
thrust  than  the  UH60  blade,  and  the  Figure  of  Merit 
is  significantly  less.  As  before,  there  is  no  flow  separa¬ 
tion  on  the  blade  except  in  the  very  tip  region.  The  tip 
vortex  location  is  slightly  inboard  of  the  UH6Q  blade 
and  has  less  tight  braiding  than  the  UH60  blade  at  a 
comparable  blade  location  although  this  blade  also  has 
a  swept-tip  planform. 

BERP  Rotor 

The  BERP  blade  has  an  unorthodox  planform 
IfiOshape  with  a  twist  distribution  along  its  length  as 
shown  in  Fig.  2.  The  details  of  the  blade  geometry 
is  described  elsewhere  (see,  for  example.  Ref.  7).  The 
main  feature  of  its  geometry  is  a  leading  edge  notch 
(see  Fig.  1)  associated  with  a  paddle-type  shape  with 
increased  chord  and  swept  delta-wing-like  planform  in 
the  tip  region.  This  unusual  planform  and  high  per¬ 
formance  airfoils  have  enabled  BERP  rotor  to  achieve 
high  performance  in  high  speed  forward  flight.  In  or¬ 
der  to  compare  planform  effects  of  this  blade  with  UH60 
blade,  two  calculations  are  made  for  collective  pitch  of 
0c  =  9®  and  13®  at  a  tip  speed  of  Mtip  =  0.628  and  Re 
=  2.95  X  10®. 

Figure  13  shows  typical  surface  pressure  for  the 
BERP  rotor  at  different  radial  locations  for  the  flow 
condition  corresponding  to  0c  =  13®.  Currently  no  ex¬ 
perimental  data  is  available  for  these  flow  conditions. 
However,  these  results  and  the  surface  pressures  for  0c 
=  9®  appear  to  be  in  good  qualitative  agreement  with 
the  nonrotating  blade  results  of  Ref.  7.  The  sectional 
thrust  variation  is  shown  for  this  blade  for  0c  =  9®  in 
Fig.  14,  along  with  UH60  results  at  identical  collec¬ 
tive  pitch  and  flow  conditions.  Since  BERP  blade  has 
variable  chord  along  its  radius  in  contrast  to  the  UH60 
blade,  the  sectional  thrust  curve  for  BERP  blade  has 
an  unusual  shape  in  the  paddle  part  of  the  blade.  If 
the  results  of  Fig.  14  are  recast  in  terms  of  equiva¬ 
lent  chord,  then  these  results  for  BERP  blade  will  have 
the  familiar  form  as  shown  in  Fig.  15.  As  seen  in  this 
comparison,  the  combination  of  leading-edge  notch  and 
increased  area  of  paddle  part  of  BERP  has  produced  an 
increased  thrust  for  the  blade  for  this  region.  Inboard, 
however,  it  produces  relatively  lesser  thrust  perhaps  be¬ 
cause  of  slightly  lower  twist  in  this  region.  For  both 


conditions  considered,  the  BERP  blade  produces  more 
thrust  at  the  expense  of  significantly  more  power. 

Table  1  summarizes  important  hover  performance 
results  for  this  rotor.  Figures  16-17  show  the  sectional 
thrust  and  bound  circulation  distributions  for  BERP 
blade  for  0c  =  9®  and  13®  at  A/„p  =  0.628  and  Re  = 
2.95  X  10®.  In  essence,  these  two  plots  show  the  same 
quantity  but  looking  in  different  ways.  The  high  thrust 
condition  produces  80%  more  thrust  at  the  expense  of 
160%  more  pov.er  compared  to  the  moderate  thrust 
condition,  corresponding  to  0c  =  9®,  at  the  same  tip 
speed.  While  the  curves  of  Figs.  15  and  16  have  the 
same  qualitative  behavior  as  that  of  UH60  blade  shown 
in  Fig.  7,  the  bound  circulation  distribution,  shown  in 
Fig.  17,  has  the  same  qualitative  behavior  as  the  sec¬ 
tional  lift  distribution  of  nonrotating  BERP  blade^.  Of 
particular  significance  is  the  distribution  in  the  vicin¬ 
ity  of  the  notch.  As  observed  for  a  nonrotating  blade 
from  both  the  wind  tunnel  tests  and  calculations^,  the 
present  calculations  also  show  that  the  sectional  thrust 
of  the  blade  in  this  region  first  decreases  and  then  in¬ 
creases  for  the  condition  of  0c  =  13®.  In  contrast,  the 
moderate  thrust  condition  results  does  not  show  this 
behavior. 

The  surface  particle  tracer  pictures  of  BERP  blade 
are  shown  in  Fig.  18  for  0c  =  9®  and  13®.  For  both 
the  thrust  conditions,  the  flow  stays  attached  on  the 
paddle  part  of  the  blade.  The  high  thrust  condition 
corresponding  to  0e  =  13®,  however,  produces  a  com¬ 
plicated  flow  structure  in  the  notch  region  identified  as 
vortical  flow  structure.  In  spite  of  the  strong  supercrit¬ 
ical  flow  for  this  condition,  as  in  UH60  blade  at  0c  = 
13®,  the  BERP  blade  produces  relatively  small  extent 
of  shock-induced  flow  separation  only  in  the  inboard  of 
notch  region  as  seen  in  Fig.  18.  The  details  of  pres¬ 
sure  contours  on  the  upper  surface  of  the  blade  and  at 
one  spanwise  location  in  this  region,  shown  in  Fig.  19, 
highlights  the  extent  of  transonic  flow  for  this  flow  con¬ 
dition.  In  contrast,  the  UH60  blade  produces  a  severe 
shock  induced  flow  separation  over  a  large  part  of  the 
blade  for  this  flow  condition.  This  total  absence  of  flow 
separation  on  the  paddle  part  of  the  blade  at  transonic 
conditions  has  a  favorable  influence  on  the  performance 
of  the  blade.  In  the  tip  region,  the  case  of  0c  =  13®  has 
produced  a  larger  extent  of  separation  on  the  upper 
surface  than  the  0c  =  9*  case.  Since  these  blades  have 
very  high  inboard  twists,  both  conditions  produce  se¬ 
vere  flow  separation.  The  extent  of  separation  for  the 
0c  =  13®  case,  shown  in  Fig.  20,  is  much  more  massive 
than  the  0c  =  9®  case  as  for  UH60  blade  case. 

Comparison  of  surface  particle  tracers  details  for 
the  IIH60  (Fig.  9)  and  BERP  (Fig.  18)  rotors  in  the 
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outer  part  of  the  blades  towards  the  tip  indicate  the 
following  features.  Since  the  blades  have  totally  dif¬ 
ferent  planforms  and  twist  distributions  in  this  region, 
they  produce  very  different  flow  patterns.  The  UH60 
blade  has  a  hook-like  twist  distribution  in  the  swept 
part  near  the  tip  while  the  BERP  blade  has  negative 
constant  twist  with  a  large  paddle  part  and  a  delta-wing 
like  tapered  tip.  For  the  low  thrust  case,  corresponding 
to  =  9°,  both  blades  show  completely  attached  flow 
except  for  the  tip  region.  In  contrast,  at  the  high  thrust 
condition  {dc  =  13“),  they  have  totally  different  surface 
flow  patterns.  The  UH60  blade  shows  strong  shock- 
induced  separation  and  reattachment  near  the  leading 
edge  area  on  _-.ther  side  of  kink  (begining  of  sweep-back 
position),  whereas  the  BERP  blade  shows  completely 
attached  flow  on  the  paddle-part  except  in  the  extreme 
tip  region.  The  notch  region  has  a  complicated  flow 
separation  and  vortical  flow  structure. 

Figure  21  shows  the  details  of  a  typical  tip  vortex 
formation  and  roll  up  process  for  the  BERP  blade.  This 
type  of  flow  visualization  picture  is  generated  by  releas¬ 
ing  unrestricted  flow  particle  tracers  in  the  tip  region  on 
both  sides  of  the  blade  at  different  chordwise  and  span- 
wise  locations  and  at  different  heights  from  the  surface 
of  the  blade.  The  particles  so  released  bunch-up  and  get 
braided  over  each  other  to  form  a  distinct  structure  of 
concentrated  vorticity.  Such  a  structure  represents  the 
formation  process  of  a  concentrated  tip  vortex  as  seen 
in  Fig.  21.  As  this  tip  vortex  is  still  forming,  it  starts 
to  roll-up  and  in  the  process  stays  distinctly  above  the 
wake  vortex  sheet.  Once  it  leaves  the  surface,  it  con- 
vects  with  the  wake.  As  this  vortex  wake  convects,  it 
contracts  in  size  and  descends  with  the  flow.  With  the 
present  nonadaptive,  structured  grids  used,  the  cap¬ 
tured  vortex  structure  is  smeared.  This  means  the  core 
is  diffuse  but  the  circulation  is  correctly  captured  in 
strength. 

Comparing  the  vortical  flow  structures  of  UH60 
and  BERP  blades  corresponding  to  0^  =  13“ ,  both  have 
distinct  but  diffuse  tip  vortex  structures.  The  notch 
area  of  BERP  produces  a  complex  vortical  stucture  as 
seen  in  Fig.  21.  Similary,  the  kink  of  UH60  blade  also 
produces  a  very  mild  vortical  structure  for  this  condi¬ 
tion.  Comparison  of  tip  vortex  formation  processes  of 
UH60  blade  and  equivalent  rectangular  blade  for  the 
low  thrust  condition,  corresponding  to  9^  =  9“  case, 
indicates  that  the  equivalent  rectangular  blade  shows 
similar  formation  and  roll-up  process  as  that  of  UH60 
blade  except  that  the  braiding  of  the  particle  tracers 
occurs  much  farther  in  the  wake  in  conformity  with  the 
observations  made  in  Ref.  15. 


Farfield  Boundary  Conditions 

A  difficult  aspect  of  numerical  prediction  of  hov¬ 
ering  rotor  flows  is  the  prescription  of  farfield  bound¬ 
ary  conditions.  In  hover,  if  the  specification  of  farfield 
boundary  conditions  is  based  on  the  quiescent  flow  out¬ 
side  of  the  computational  box,  then  the  flow  into  and 
out  of  the  box  will  be  zero.  This  creates  a  closed  box 
atmosphere  for  the  rotor,  with  the  flow  inside  the  box 
recirculating  when  the  rotor  is  spinning.  In  princi¬ 
ple,  this  is  precisely  the  environment  in  many  hover 
test  chambers.  However,  a  large  computational  box 
of  this  size  is  not  cost  effective  and  exceeds  the  mem¬ 
ory  limits  of  present  day  supercomputers.  Alterna¬ 
tively,  with  a  smaller  computational  box  around  the 
rotor  one  should  specify  non-zero  flow  farfield  bound¬ 
ary  conditions  that  will  allow  the  flow  to  enter  and 
exit  the  box  without  violating  the  conservation  laws. 
Currently,  no  precise  methods  are  available  to  spec¬ 
ify  these  boundary  conditions,  although  some  schemes 
have  been  suggested^.  The  present  calculations  take  a 
slighty  different  approach  to  specifying  farfield  condi¬ 
tions.  A  three-dimensional  sink  whose  strength  is  de¬ 
termined  from  the  thrust  of  the  rotor,  is  located  at  the 
axis  of  rotation  in  the  plane  of  the  rotor.  With  this,  the 
flow  enters  the  computational  box  from  all  boundaries, 
except  through  an  exit  plane  on  the  lower  boundary 
whose  area  is  half  the  area  of  the  rotor  disk.  Across 
this  exit  plane,  a  characteristic-type  outflow  boundary 
condition  is  specified  by  prescribing  the  w-velocity  as 
twice  its  value  at  the  rotor  plane  from  the  momentum 
theory  concept.  The  other  four  quantities  are  extrapo¬ 
lated  from  within. 

Using  the  above  approach  of  specifying  the  farfield 
boundary  conditions,  calculations  were  repeated  for  the 
UH60  blade  at  flow  conditions  of  Mup  =  0.628,  9c  =  9“ 
and  Re  =  2.75x10®.  Figure  22  shows  a  typical  result  of 
sectional  thrust  distribution  compared  with  that  gener¬ 
ated  with  old  boundary  conditions  precedure  and  with 
experiments.  Although  there  is  general  agreement  of 
overall  results  in  this  plot,  the  reasons  for  small  dif¬ 
ferences  seen  near  the  kink  and  in  the  tip  region  is 
not  clear  at  this  time.  The  performance  parameters 
are  shown  in  Table  1.  As  seen  in  this,  the  thrust  and 
power  values  are  about  3%  more  than  that  of  the  old 
boundary  conditions. 

The  important  test  of  this  new  boundary  condition 
study  is  the  behavior  of  the  captured  wake  trajectory. 
After  identifying  the  fluid  particle  tracers  that  reside 
with  the  core  of  the  tip  vortex,  such  as  shown  in  Fig. 
21,  one  or  more  fluid  particle  tracers  were  released  in 
the  tip  region  of  each  blade  and  the  evolution  of  these 
tracers  in  time  was  monitored.  Figure  23  shows  two 


6 


views,  a)  contraction,  and  b)  descent,  of  such  a  wake 
trajectory  for  the  above  test  case.  The  nearfield  tra¬ 
jectory  of  this  wake,  for  both  contraction  and  descent, 
appears  to  be  in  reasonable  agreement  with  the  exper¬ 
iments  and  is  shown  in  Fig.  24.  Previous  calculations, 
for  the  same  computational  box  with  old  boundary  con¬ 
ditions,  showed  that  the  wake  trajectory  for  this  four- 
bladed  rotor  was  accurately  captured  only  up  to  about 
180®  of  vortex  age  (azimuth  travel).  Subsequently,  the 
wake  trajectory  continued  to  descend  but  expanded  in 
size,  from  its  previous  size  at  180®  azimuth,  and  eventu¬ 
ally  got  cought  in  the  recirculating  flow.  With  the  new 
farfield  boundary  conditions,  the  flow  exits  smoothly  at 
the  bottom  boundary  and  so  does  the  wake. 

Concluding  Remarks 

A  limited  airloads  and  flowfield  analysis  of  the 
UH60  model  rotor  and  the  BERP  rotor  is  presented 
for  lifting  hover  configurations.  The  calculation  pro¬ 
cedure  uses  a  computational  fluid  dynamics  Navier- 
Stokes  capability  to  capture  the  flowfield  and  nearfield 
vortical  wake  and  its  trajectory  without  having  to  pre¬ 
scribe  any  wake  models.  At  present,  the  vortical  struc¬ 
ture  in  the  wake  is  smeared  due  to  grid  coarseness,  and 
this  aspect  of  the  method  needs  further  improvements. 
Nevertheless,  the  induced  flow  in  the  plane  of  the  rotor 
and  the  near-field  wake  trajectory  are  computed  fairly 
accurately,  and  the  calculated  pressures  on  the  blade 
agree  well  with  experiments. 

The  results  and  discussion  presented  here  focus  on 
evaluating  the  importance  and  influence  of  planform  ef¬ 
fects  on  the  airloads  and  vortex  wake  geometry.  The 
results  for  the  UH60  rotor  show  good  agreement  with 
experimental  data  for  the  surface  pressures,  thrust, 
power  and  fair  agreement  for  wake  trajectory.  For  a 
given  flow  condition,  the  BERP  rotor  produces  more 
thrust  at  the  expense  of  increased  power  and  reduced 
Figure  of  Merit.  The  results  for  the  UH60  and  BERP 
blades  at  transonic  high  thrust  conditions  indicate  that 
the  BERP  blade  produces  milder  shock-induced  sepa¬ 
ration  than  the  UH60  blade.  The  BERP  blade  also 
shows  tightly  braided  tip  vortex  structure  compared  to 
the  UH60  blade.  The  good  comparison  of  the  UH60 
rotor  performace  calculations  with  experiments  indi¬ 
cates  that  this  free-wake  Navier-Stokes  CFD  capability 
is  promising  as  a  tool  for  analysing  the  properties  of 
new  and  exotic  blade  shapes  whose  properties  are  not 
known  a  priori. 
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Table  1 

Calculated  and  Experimental  Hover  Performance  Parameters  for 
UH-60  Blackhawk  and  BERP  Rotors  for  M,.  =  0.628 


Rotor 

Methods  of  results 

0c 

Cyo 

Cq/o 

FM 

UH-60  model  rotor 

Experiments  (Ref.  6) 

Unknown 

■a 

0.0069 

0.73 

N-S  results 

9° 

0.084 

0.0068 

0.73 

N-S  results  with  new  farfield  BC 

9° 

0.087 

0.0069 

0.76 

UH-60  equivalent 
rectangular  blade  rotor 

N-S  results 

■ 

0.078 

0.0065 

0.68 

UH-60  linear 
twist  biade  rotor 

N-S  results 

■ 

m 

0.0072 

m 

UH-60  model  rotor 

N-S  results 

13° 

0.140 

0.0149 

0.72 

BERP  rotor 

N-S  results 

9° 

0.094 

0.0098 

0.68 

N-S  results 

13° 

0.168 

0.0259 

0.62 
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Fig.  5  Surface  pressure  distributions  at  various  span  stations  for  the  UH60  blade  compared  with  experi¬ 
mental  data.  Mtip  —  0.628,  Re  =  2.75  x  10®,  and  calculated  CtI<^  =  0.084. 


Fig.  6  Sectional  thrust  loading  for  the  UH60  blade, 
equivalent  UH60  rectangular  blade,  and  the  UH60 
blade  with  linear  twist  distribution  in  the  tip  region. 
Mtip  =  0.628,  =  9®,  and  Re  =  2.75  x  10®. 
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Fig.  7  Comparison  of  sectional  thrust  loading  for  the 
UH60  blade  for  two  thrust  conditions.  Mup  =  0.628, 
and  Re  =  2.75  x  10®. 
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Fig.  8  Upper  surface  particle  flow  (skin  friction)  patterns  of  the  UH60  blade  for  Be  =  9“  and  13®  cases  at 
Mup  =  0.628  at  Re  =  2.75  .x  10®. 
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Fig.  9  Close-up  view  of  surface  particle  flow  in  the  tip  region  for  the  UH60  blade  for  6c  =  9®  and 
9c  =  13®  at  Mtip  =  0.628,  and  Re  =  2.75  x  10®. 
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Fig.  10  Mach  contours  for  the  UH60  blade  for  the  conditions  of  Mup  =  0.628,  Oc  =  13®,  and  Re  =  2.75 
X  10®.  Views  show  a)  upper  surface  contours  at  L  =  8  (away  from  the  wall),  and  b)  cross-sectional  plane 
at  r/R  =  0.94. 


Fig.  11  Comparison  of  calculated  and  experimental  Fig.  12  Comparison  of  calculated  and  measured  Figure 

hover  performance  for  the  UH60  model  rotor.  of  Merit  for  the  UH60  model  rotor. 
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Fig.  13  Typical  surface  pressure  distributions  for  the  BERP  rotor  blade  at  different  span  stations;  Mup 
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Fig.  14  Radial  thrust  loadings  (based  on  local  chord) 
for  UH60  and  BERP  blades  at  the  same  collective  pitch 
setting  of  0c  =  9®  with  Mup  =  0.628. 
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Fig.  15  Radial  thrust  loadings  (based  on  equivalent 
chord)  for  UH60  and  BERP  blades  at  the  same  collec¬ 
tive  pitch  setting  of  9e  =  9®  with  Mup  =  0.628. 
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Fig.  16  Comparison  of  radial  thrust  loadings  for  BERP 
blade  at  =  9“  and  Oc  =  13“  with  M^p  =  0.628  and 
Re  =  2.94  X  10®. 


Fig.  17  Bound  circulation  distributions  for  BERP  blade 
for  the  conditions  of  Fig.  16. 
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Fig.  18  Surface  particle  flow  details  mimicing  oil  flow  patterns  for  the  BERP  blade  at  6c  =  9®  and  13“; 
Mtip  =  0.628  and  Re  =  2.94  x  10®. 
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b)  View  at  section  A-A 


a)  Upper  surface  pressure  contours 


Fig.  19  Pressure  contours  of  the  BERP  blade  at  Mup  =  0.628,  9c  =  13®,  and  Re  =  2.95x10®.  Views  show 
a)  upper  surface  and  b)  cross-sectional  plane  at  r/R  =  0.80. 
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Fig.  20  Upper  surface  particle  flow  tracers  pattern  for  the  BERP  blade  at  Mup  =  0.628,  9c  =  13®,  and 
Re  =  2.94  X  10®. 
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Tip  vortex  formation 

Fig.  21  Near-field  view  of  tip  vortex  formation  ana 
roll-up  process  for  the  BERP  blade  at  6e  =  13“,  Mup 
=  0.628. 


Fig.  22  Sectional  thrust  distributions  for  the  UH60 
blade  with  old  and  new  farfield  boundary  conditions 
compared  with  experimental  data.  Mtip  =  0.628,  9c  = 
9®  and  Re  =  2.75  x  10®. 
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(a)  Contraction  (b)  Descent 

Fig.  23  Calculated  tip  vortex  trajectory  for  a  4-bladed  UH60  model  rotor  in  hover  showing  wake  contrac¬ 
tion  and  descent;  Mup  =  0.628,  9c  =  9.0®,  Re  =  2.75  x  10®. 


Fig.  24  Comparison  of  wake  contraction  and  descent  for  a  4-bladed  UH60  model  rotor  with  experimental 
data  for  the  conditions  of  Fig.  22. 


